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Abstract. We consider the asymmetric random average process (ARAP) with 
continuous mass variables and parallel discrete time dynamics studied recently by 
Krug/Garcia and Rajesh/Majumdar [both Jrl. Stat. Phys. 99 (2000)]. The model is 
defined by an arbitrary state-independent fraction density function 4>(f) with sup- 
port on the unit interval. We examine the exactness of mean-field steady-state mass 
distributions in dependence of <j) and identify as a conjecture based on high order 
calculations the class M of density functions yielding product measure solutions. 
Additionally the exact form of the associated mass distributions P(m) is derived. 
Using these results we show examplary the exactness of the mean-field ansatz for 
monomial fraction densities 4>(r) = (n — l)r' 1-2 with n > 2. For verification wc 
calculate the mass distributions P(m) explicitly and prove directly that product 
measure holds. Furthermore we show that even in cases where the steady state is 
not given by a product measure very accurate approximants can be found in the 
class M. 

Keywords: Non-equilibrium physics, stochastic systems, interacting particle sys- 
tem, random average process, q model, invariant product measure, discrete time 
dynamics, exact solution 



1. Introduction 

Interacting particle systems far from equilibrium represent due to their 
wide applications in physics and other related topics like traffic flow 
modeling [1] or force propagation in granular media [5] a popular the- 
oretical research field [8]. In this paper we focus on a model studied 
recently by Krug and Garcia [4] and Rajesh and Majumdar [7] that 
is closely related to the q model introduced by Coppersmith et al. 
[3]. It describes the movement of particles on the real line according 
to a given probability distribution depending on the distance of two 
particles. This system is equivalent to a stick model where the height 
of the sticks represents the particle gap [4, 7]. Below we present a brief 
description for completeness and demarcation of our framework. 

The asymmetric random average process (ARAP) is defined on a 
one-dimensional periodic lattice with L sites. Each site i carries a non- 
negative continuous mass variable rrij. In every time step t — > t + 1 
for each site a random number rj is generated from a time- and site- 
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independent probability distribution 4> defined on [0,1]. The fraction 
Ti determines the amount of mass r^mj transported from site i to site 
i + 1 (asymmetric shift). 

In the following we concentrate on the infinite system in the thermo- 
dynamic limit, i.e. L — > oo and M = J2i m i ~~ * 00 with finite constant 
density p = ^. Furthermore we examine only time independent steady 
state dynamics. 

In case of a uniform fraction density <j){r) = 1 mass distributions 
and moments have been calculated in mean-field approximation for 
different types of dynamics [4, 6, 7]. Although these analytical results 
show excellent agreement with numerical simulations, one can prove 
that for random sequential update a product measure ansatz fails [7]. 
For the fully parallel update, however, one can adopt the result of [3] 
where the exactness of the mean-field approach is proven in the context 
of force fluctuations in bead packs described by the q model [3, 2]. 

In this paper we like to reinforce investigations by considering ar- 
bitrary state independent fraction density functions (j) = <p(r) and 
focussing on the fully parallel update only. Such density functions 
are relevant for practical applications, e.g. as approximation for the 
ARAP with a finite mass cut-off [11] or for describing suitable bead 
pack problems. 

In section 2 a functional equation acting on the Laplace-space of the 
single-site mass distributions P(m) is derived. Product measure holds 
if a mean-field ansatz is a solution of this condition. 

In section 3 we determine the set M of all fraction probability 
densities <p(r) yielding product measure states. This result represents 
a conjecture based on exact high order calculations. Now statements 
about the exactness of mean-field are possible without solving the 
master equation. We also derive the corresponding mass distributions 
P(m). 

As an example the ARAP with monomial ^-function is studied (sec- 
tion 4). We calculate P(m) and prove the exactness of mean-field by 
using the explicit form of P(m) on the one hand and the criterion of the 
last section on the other hand. Due to the fact that the n dimensional q 
model with uniform distributed c/'s leads to the same mass distribution 
[3] its relationship with the ARAP is discussed, too. 

In section 5 we demonstrate that the results obtained in the previous 
sections can be used to obtain very accurate approximations for the 
mass distributions of ARAPs which do not lead to product measures. 

In section 6 we conclude with a summary and give a brief outlook 
for further work. 
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2. A Functional Equation of Exact Mean-Field Solutions 



In this section we derive a functional equation (see eq. (7) below) for 
determining and testing mean-field solutions. 

The fundamental element of all upcoming considerations is the mas- 
ter equation. For the single site mass distribution P(m) in the station- 
ary state it reads [7] 

/•oo /*oo rl rl 

P(m 2 ) = / dmi I dm,2 P(mi,m2) / dr\ I dr 2 <Kn)<K?"2) 
Jo Jo Jo Jo 

x 5 (m' 2 - [rpi + (1 - r 2 )m 2 ]) . (1) 

The (5-function ensures mass conservation. We have assumed transla- 
tional invariance so that the distribution is site- independent. 

By the mean-field ansatz P(mi, mi, ■ ■ ■) = Y\iP(fni) equation (1) 
determines the single site distribution P(m). The resulting product 
measure is exact if the mean-field ansatz holds for all joint probabilities 
in the stationary state, too, i.e. P(m) has to satisfy 




x J| S (m- - [n-mii-! + (1 - ri)mi}) (2) 

for all k e 1N>2- For k = 2 this equation reduces to (1). This ap- 
pears to be an infinite set of conditions, but Laplace transforming 
P(mi, mi-, ■ ■ •) reduces (2) to just one functional equation. By intro- 
ducing the k— dimensional Laplace-transform 

/•oo 

Q( Sl ,...,s k )= d k mP(m)e~^ s \ (3) 
Jo 

where (m, s) = J2i=i miSi, and using the map 

F Q (s, §) = Cdr 4>{r)Q ((1 - r) s + rs) (4) 
Jo 

equation (2) reads in Laplace space 

k k-1 

HQ ( Si ) = F Q (0, Sl ) ■ H F Q ( Si ,s i+1 ) ■ F Q (s k ,0) (5) 

i=l i=l 

for k e iV. 
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By a straightforward proof we now show that the conditions for 
k 7^ 2 are redundant. The k = 1-equation 

Q(s) = F Q (s,0)-F Q (0,s) (6) 

is used to determine Q. We rewrite the k = 2-criterion using (6) and 
obtain 

Fq(s 1 ,S2) = Fq(s 1 ,0)-F q (0,s 2 ) . (7) 

Applying (6) and (7) proves the validity of (5) for all k > 3 and using 
the identity Q(s) = Fq(s, s) equation (6) becomes a special case of (7). 
So (7) is the only necessary equation to determine a mean-field solution 
(si = S2) and check its accuracy (s\ / 52)- 



3. Exact Mean-Field Solutions 

In this section we derive the set of density functions (p(r) that result in 
product measure steady states. This yields a more useful criterion for 
determining the exactness of a mean-field solution without calculating 
and verifying Q(s) by condition (7). In addition we derive the mass 
distribution P(m). 

We start by proving that equation (6) has always a unique solution 
in the space of functions that are analytical in the origin. Intuitively 
one supposes this feature because the mass moments 

m n = {m n ) p{m) = / dmm n P(m) (8) 
J 

are (formally) generated by m n = (— l) n Q^ n \o). But a priori we do 
not know if all derivatives of the moment function Q(s) exist or 
ensure convergence. The moments mo = 1 and mi = p are determined 
by the normalization and the density p, respectively. 

We first represent the moment function as a (formal) power series, 
i.e. Q(s) = J2 n a nS n . The coefficients a n are related to the moments (8) 
by a n = (— l) n! ^f and thus we have do = 1 and ai = —p. The remaining 
coefficients are determined with help of a recurrence relation obtained 
by inserting the series representation into (6): 

n— 1 

a n = z Vk,o Po,n-k a k a n ~k (Vn > 2). (9) 

1 - Pn,0 ~ P0,n k=1 

Here p n ,m are generalized moments of the fraction density cp defined by 
Pn, m = (r n (l - r) m ) Hr) = jf 1 dr<f>(r)r n (l - r) m . (10) 
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We assume /x n ,m > in the following which is equivalent to 4>{r) / 
for r £ (0, 1). ix n ,m = does not occur for continuous distributions, but 
e.g. for (p(r) = p5{r) + (1 — p)5(l — r). This ARAP is trivial for p = 
or p = 1 and not solvable under mean- field assumptions for < p < 1. 
From 

(10) n_1 fn\ 

1 - Mn,0 - A*0,n = Xj \ ky k,n ~ k > (Vra > 2) (11) 

we conclude that all a n are well-defined and the solution of (6) is 
unique. By the formula of Cauchy-Hadamard we then show that Q 
is holomorphic in s = 0: We start by proving the lemma 

Kl < D^cnp 71 (Vn > 1) (12) 

with D = 1 _ ^_ - and the density p. Here {cn} nG ^y are ^ ne Catalan 
numbers [10] fulfilling the equations 

ci = l and c n = ^c fc c n _ fc = - l " M (Vn > 2). (13) 

Inserting (12) into (9) using (13) and the fact that /x n>m > H n +k 1 ,m+k 2 
for all fcj 6 IN shows inductively the validity of the lemma (12). From 
linin^oo ^fc^ = 4 we conclude that the series expansion Q(s) = J2 n a nS n 
has a positive radius of convergence. 

After proving that mean-field solutions are always representable as 
power series we now try to express the exact solutions in terms of the 
density function <f>. Inserting Q(s) = J2n a nS n into (7) yields an infinite 
set of conditions 

{^j Vk,n-k a n = Wfc,0 fM),n-k «fe «n-fe (Vn G Wo, V/c G {0, . . . , n}). 

(14) 

In the following we refer to a special equation of (14) as (n, /c)-condition 
or equation of order n. Summing over k = 1, . . . , n— 1 in (14) yields con- 
dition (9). This shows again that (7) includes formula (6) and implies 
convergence of Q(s). 

Now we are confronted with the problem that a n is determined by 
n + 1 equations. Since the (n, 0)- and (n, n)-conditions match identities 
there are effectively n — 1 equations to be fulfilled for n > 2. Thus for 
n > 3 the occurence of inconsistencies is possible and for arbitrary <j) 
or moments n n ,m we see by explicit calculation contradictions in order 
n = 3 already. Is it possible to find a set of moments {/J, n<m } such that 
(14) is satisfied for all (n, fc)? 
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Assuming that all (n, fc)-conditions yield the same a n we see from 
(14) that the function 

f(n, k) = 7HT- • (15) 

\k) Hk,n-k 

is independent of k. Therefore the n — 2 consistency equations 

f(n, 1) = f(n, 2) = ... = f(n, n-l) = a n (16) 
have to be satisfied. The definition (10) implies that fJ, n ,m can 

be 

expressed by fijfi only: 



Hn,m = Y^[ ■ \(- 1 ) 3 with Hj = Hj,0- (17) 




Together with (14) it follows inductively that a n = a n ([Xo, . . . , fi n ). 
This leads to a successive constructive approach in n: We try to find 
a fi n (depending on fiQ, . . . , Hn-i) that solves all conditions (16) of 
order n starting with n = 3 and repeat this for all upcoming orders 
n = 4,5,... . 

The n — 2 equations (16) are linear in (i n . This ensures uniqueness 
of a possible solution. For arbitrary k, k = 1, . . . ,n — 1 with k / k we 
obtain 

_ (fc) ^ n < k 9 n ,k ~ (fc) ^n,k 9n,k 

Here g ntk = fJ-kfi Mo,n-fc afc a«-fc and = E™=o fc ) ("l^ Mfc+j 

only depend on /xi, . . . , fin-i- Thus (18) gives us the desired recur- 
sion relation to determine all moments \x n . However, we have to check 
whether the r.h.s. of (18) is independent of k and k. We have done this 
up to order n = 10 by computer algebra and conjecture the validity for 
all n. Our results furthermore lead us to conjecture the following form 
of the solution of the density moments: 

_ n TT 1 Z + A 1 _ T(n + X 1 ) T(A 2 ) 

tk l + >*~ r(Ai) r(n + A 2 ) liyj 

with 

A l — Ml o ' A 2 — 2 ' \ Z(] > 

H2 -HI Hi- Ml 

which again has been checked up to n = 10. Note that jjli and ^ 
are free parameters that only have to be choosen with respect to the 
general moment properties 

1 > hi > M2 > Mi ■ (21) 
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The special case p,i = /Lt 2 yields fj, n = p r { representing <p(r) = S(i — fi±) 
which leads to Q(s) = exp(— ps) and thus to the mass distribution 
P(m) = 5(m — p). So we assume /i 2 > Mi i n the following. Under these 
restrictions p, n >3 is pole-free and satisfies < p n +i < Mn as demanded. 

(19)-(21) define the set M of all fraction densities <j)(r) yielding prod- 
uct measure steady-state distributions. So equations (19-20) represent a 
powerful criterion for determining the exactness of a mean-field ansatz: 
We only have to calculate the moments p n of the fraction density 
(if they are not already given) and check consistency with (19) and 
(20) - without even calculating the mean-field mass distribution or its 
Laplace-transform! 

Note that the parametrization of M. in terms of p\ and /i 2 is arbi- 
trary and a consequence of our construction - other parametrizations 
are possible. For symmetric densities <j>{r) = 0(1 — r) the space A4 
reduces to one dimension because \i\ = \ is fixed. 

After determining the class M of fraction densities leading to a prod- 
uct measure we now like to calculate the single site mass distribution 
P(m) for these <f> G M. From the (n + 1, n)-condition (14) we derive 
the recurrence relation 

A 2 (n + lK +1 +p(n + A 2 K = (22) 

that is valid for all n G JNq. (22) corresponds to a first order differential 
equation for the moment function Q(s) = J2 n a nS n - 

(A 2 + ps) Q'(s) + \ 2 p Q{s) = . (23) 

Using the boundary condition Q(0) = 1 we obtain 



Q(s) = 7 — 7v^ (24) 



or, by calculating the inverse Laplace-transform, 

In contrast to fi n depending on both Ai and A 2 - see equation (19) - 
P(m) is a function of A 2 only. So ARAPs with A 2 fixed and Ai arbitrary 
have identical mass distributions (25). 
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4. Explicit solution for monomial density function <j>(r) 

In this section we derive the solution of the ARAP with density function 

Mr) = (n- l)r n - 2 (Vn G M> 2 ) (26) 

in a closed form (n — 1 is the normalization constant) and prove the 
exactness of the product measure Yik P{ m k) both explicitly and using 
the criterion (19-20). Additionally a brief comparison between ARAP 
and the q model is presented. 

We start by constructing the analytic solution of the functional 
equation (6) explicitly. For n = 2, where (f> n reduces to a constant 
distribution, we refer to [3, 4, 7] and find 

w =otw (27) 

To solve the n = 3-problem we generalize the method used in [3, 4, 7]. 
Defining the functions 

V(s) = f 1 drQ 3 (rs) and W(s) = C dr r Q 3 (rs) (28) 
Jo Jo 

the functional equation (6) transforms into 

Q 3 (s) = 4W(s) (V(s) - W(s)) . (29) 
From (28) we derive 

sV'(s) + V(s) = Q 3 (s) = sW'(s) + 2W(s) (30) 
which implies the following relation between V and W: 

V(s) = W{s) + i J W{s)ds . (31) 

Defining the antiderivative f(s) = fW(s)ds and inserting (31) and 
(30) into (29) yields the nonlinear differential equation 

s 2 f"(s) + 2sf( S )-4f'(s)f( S ) = (32) 

with boundary conditions /'(0) = \ and /"(0) = — This results in 
f(s) = | (1 + |s) 1 , e.g. solved by power series ansatz, and we get 
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The results (27) and (33) suggest the assumption 

= (i + E s f ' < 34 ) 

for general n. Q n fulfills the initial conditions Q n (0) = 1 and Q' n (0) = 
—p. By a straightforward induction in n (using partial integration) we 
are able to prove 

%(^) = ttp; and F Q"(M= * . w -i ( 35 ) 

i+ n S i 1 +n S ) 

and see with (34) that (6) is valid. 

The next step is to verify the functional equation (7). This is done 
again straightforwardly by induction in n. So (34) represents the exact 
solution of the ARAP with fraction density (26). 

Now we rederive (34) using the mean-field criterion (19-20). Calcu- 
lating the moments of (f> n exactly is an easy task and yields 

= n -i 1 , • ( 36 ) 
n — 1 + k 

Especially we have [i\ = rk ^- and fi2 = ^qrj 1 and, using (20), we obtain 
Ai = n — 1 and A2 = n. 

The exact form (36) of the moments is reproduced by taking into 
account (19). Thus the monomial density functions <p n are elements of 
A4. Using A2 = n shows the equivalence of (24) and (34) after all. 

Additionally this example confirms the validity of the conjecture 
(19-20). 

For completeness we also give the explicit form of the single-site 
mass distribution for the monomial density functions (26): 

Finally we like to mention that Q n also satisfies the relation 



Qn(s) = (F Qn (0,s)Y 



[ dr<p(r)Q n (rs) 
Jo 



(38) 



This functional equation represents the master equation of the n an- 
cestor q model with uniform distributed g's and was explicitly solved in 
[3]. Because of the formal difference of the underlying equations (6) and 
(38) the coincidence of the corresponding solution (34) is remarkable. 
Otherwise this result is pointing to a deeper relationship between n 
dimensional q model and ARAP. 
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As already noticed in [4] the two ancestors q model (n = 2) corre- 
sponds to an ARAP with symmetric density function 4>{r) = 0(1 — r). 
Coppersmith et al. [3] identified a set of mean-field solutions gener- 
ated by monomial distributions, i.e. f(q%j) = q m , yielding <j){r) = 
^ra+s) ! r m (l — r) m and it is easy to see that these densities are also 
elements of M. 



5. Approximative mass distributions for arbitrary density 

functions 

As an application of our results, we construct approximative mass dis- 
tributions for arbitrary density functions <p. This is done by calculating 
the parameter A2 = A2(</>) with the help of (20) (using the exact mo- 
ments Hi and hi of 4>) and taking the corresponding mean-field solution 
(25) as an approximation. To illustrate this method and estimate its 
quality we consider two examples. 

Our discussion starts with density functions being convex combina- 
tions of elements of M . We restrict here to the special case 

<\> c = (1 - c)02 + c^3 , c £ [0, 1] (39) 

with the monomial density functions <p n defined in (26). This convex 
combination of probability densities conserves their basic properties 
like normalization or positivity. Calculating the first and second mo- 
ment of (f) c yields A2 = • Inserting this result into (25) generates 
c-dependent approximations P c . 

Comparing the distribution P c with numerical data shows an excel- 
lent agreement between approximation and the results of Monte Carlo 
simulations for all values of c (see figure 1). Only for small masses m 
systematical differences occur. 

Furthermore one can prove that </> c M. for all < c < 1. This is 
most easily seen by comparing the third moment ^3 of (39) with the 
corresponding mean- field expression (19). So the excellent agreement 
of the data match is far from trivial. Nevertheless 4> c is an interpola- 
tion between exact mean-field solutions and may inherit some of their 
properties. 

Our second example represents the simplest version of an ARAP 
with cut-off. It is based on the model with uniform fraction density, 
but enhanced by an additional parameter r$ 6 [0,1]. The cut-off ro 
controls the movement in the following way: A stick fragment rirrii is 
only transferred if r« < ro- The corresponding density function takes 
the form: 

<f> ro (r) = (l-r )5(r) + H(r -r) (40) 
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where H(r) is the Heaviside step-function. A detailed discussion of 
several truncated ARAPs can be found in [11]. 

If 7*o = 1 we obtain the free ARAP which is exactly solvable by 
product measure ansatz [3] . For < ro < 1 one can show that (p ro M. . 
In the case ro = 0, where no motion is possible, the mean-field condition 
(7) reduces to an identity - so any distribution represents an exact 
mean-field solution. In this sense (40) interpolates between ARAPs 
with product measure steady state as in the first example, but the 
construction is not a convex combination. 

The approximants P ro are calculated as described above and match 
the Monte-Carlo simulation data perfectly again except deviations for 
small m (see figure 1). 

6. Conclusions 

We have studied the asymmetric random average process with arbitrary 
fraction density (f>. Based on the analysis of the functional equation (7) 
we have identified all ARAPs where the stationary mass distribution 
P(m) is given by a product measure. The corresponding (^-functions 
form the mean-field class Ai. They are given by their moments fi n = 
(r ra )0 which depend on two free parameters, Ai and A2, related to the 
first moments \i\ and (12 (20). These can be chosen arbitrarily subject 
only to the general conditions (21). For the product measure ARAPs 
also the stationary mass distribution P(m) can be calculated explicitly 
(25). Surprisingly it depends only on the parameter A2. 

The presented approach is rigorous except for a formal proof of (19) 
which we have conjectured on the basis of computer algebra calcula- 
tions. 

As shown our results can be used to obtain accurate approximations 
for ARAPs which do not lead to exact product measures: instead of 
calculating mean-field approximations by (6) one can use the corre- 
sponding distribution (25) with appropriately determined A2 = A2 (<?!>) 
as a first guess. A detailed discussion of M could explain why this 
method often matches perfectly with numerical simulation data. We 
suppose that arbitrary choosen <f> lie "near" to the subspace M. with 
respect to a suitable chosen norm. 

We like to mention that our approach can easily be adopted to 
the ARAP with random sequential update. Since a mean-field ansatz 
breaks down even for the simple model with uniform (^-function [7] we 
have not done calculations for this case. Nevertheless the existence 
of product measure ARAPs with continuous time dynamics is also 
possible. Furthermore one could try to transfer our calculations to the 
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q model to supplement the findings of [3] . It would be interesting to use 
our results for the calculation of other properties, too, e.g. the tracer 
diffusion coefficient [9]. 

The calculation of exact solutions of ARAPs with arbitrary <p(r) 
seems to be the final step of work. Vanishing of the 2-site-correlation 
functions [4] and the excellent quality of mean-field approximations 
have already pointed to a distribution that deviates only slightly from 
product measure form. Supplementary our work gives by equation (25) 
a rough shape of the general solutions. 

Finally it would be interesting to extend our considerations to state 
dependent density functions, i.e. <f> = (f>(r, {mj}), and especially to 
local probabilities 4>(r,m). These models are for example interesting 
for truncated asymmetric random average processes [11]. 
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Figure captions (Zielen, Schadschneider, JSP 101-110) 



Figure 1. Analytical (— ) and numerical (0, °) mass distributions 
P(m) of the convex combined ARAP (c = 0.5, left diagram) and the 
truncated ARAP (tq = 0.5 (0) and tq = 0.1 (o), right diagram). The 
analytical curves are appropriate approximants taken from the mean- 
field class M.. The numerical results are obtained by Monte-Carlo 
simulations of systems with size L = 1000 and random initial condition. 
After 10 4 steps the distribution was measured for 10 7 timesteps. A log- 
log plot is used to exhibit the deviations for small masses m which can 
hardly be seen in a conventional representation. 
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Figure 1. 
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